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^^ ■ To study the time decay laws (tdl) of quasibounded hamiltonian systems we 

f— ^ ■ have considered two finite potential wells with oscillating walls filled by non 
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Q . characteristic of trapped periodic orbits. However, the quantum dynamics do 
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I. INTRODUCTION 

The time decay laws (tdl) of classical quasibounded hamiltonian systems use to include 
meaningful information related to the bounded transient dynamics. The feature of the tdl 
corresponding to a fully chaotic system has been well established. In a special way, the 
tdl's are greatly sensible to the characteristic of the trapped periodic orbits. Systems having 
a fully hyperbolic set of trapped periodic orbits show a purely exponential tdl, while the 
tdl becomes algebraic if they have a parabolic subset of them. So, in the case of the Sinai 
billiard []T| and the Bunimovich stadium [Q , paradigms of fully chaotic hamiltonian systems, 
the existence of the so called bouncing ball orbits (parabolic and non isolated periodic orbits) 
gives algebraic long time tails in the tdl's ^^■ 

On the other hand, it is a well known fact that the characteristic and the distribution 
of periodic orbits in bounded systems are relevant subjects to quantify semiclassically such 
systems [||. 

To establish the influence of the periodic orbits in the decay of quantum systems, in 
the present work we have studied the decay processes in two quasibounded systems whose 
trapped periodic orbits have different characteristics. 

We show that the classical systems having trapped periodic orbits of different charac- 
teristics have qualitatively different tdl, while the quantum analogues do not show such 
sensitivity. In other words, we show that two classical systems with exponential and alge- 
braic tdl respectively they have algebraic quantum tdl (qtdl) when they are considered as 
quantum systems. 

Our work is organized in the following way. In Sec. || we introduce the classical systems 
whose tdl we study in Sec. |T|. Sec. |^ is devoted to describe the quantum systems and 



the numerical method that we have employed to compute the qtdl. In Sec. we show the 



resulting qtdl. Finally, Sec. VI is devoted to conclusions. We include one appendix to show 



some characteristics of the tdl's. 



II. THE CLASSICAL SYSTEMS 

In the present work we deal with a kind of one- dimensional periodically driven hamil- 
tonians. The corresponding system can be described in the extended phase space as a 
conservative system with two degrees of freedom. 

Let us consider non interacting point particles of unity mass moving inside a one- 
dimensional time dependent potential well 

oo if g < , 
Viq,t) = { ifO<g<a(t) , (2.1) 

oo if a(t) < q . 

That is a perfectly reflecting one-dimensional box with an oscillating wall according to 
the time dependent law a(t). 

Looking at the graph of the position of one particle q, as a function of time t, we can see 
that the system is equivalent to a particle forwards in a two dimensional infinite pipe with 
appropriate reflecting conditions at the boundary (See Fig. |I](a)) . The reflecting condition 
(the change of sign of the relative velocity between the walls and the particle) imposes: 

Vf = —Vi at o = , 

^ (2.2) 

Vf = 2d(t) — Vi at g = a(t) , 

where Vf{vi) is the velocity of the particle after (before) the collision. Both velocities 
correspond to the sloops of the straight lines in the q.vs.t graphs [v = tana). 

Let be a(t) periodic in time, so a{t) = a{t + r). 

Taking account the time lattice translation symmetry we can consider only an elemental 
segment of the channel setting periodic boundary conditions, namely t = r— ^t = Oso 
that q{t = r) ^ q{t = 0) and vit = t) ^ v(t = 0). So, the system can be seen as a 
two dimensional billiard where particles leaving at t = r emerge at t = having the same 
velocity v and the same position q while the collisions at the walls follows the laws (2^). In 



this version of the system, we have two degrees of freedom associated to the coordinates q 



and 9 = t — [t/T]T (reduced time) where [...] means the integer part (See Fig. |T] (b)). The 
conjugate momenta are v and —E (kinetic energy) respectively. 

When the potential well is given by rule ( p.l| ) particles inside the box are bounded despite 
their velocities. On the other hand, if the oscillating wall (at g = a(t)) is a finite barrier 



of potential of height Vq, particles with velocities | f |> V^Vq will leave the well when 
they reach the wall. In this case, according to the characteristic of the dynamics, which is 
controlled by the law a{t) of the oscillating wall, the system can be transiently bounded. 

We have focused our attention in two particular rules a{t). In the following we take r as 
a unit of time and the mean value of a{t), < a >t= 1/t Jq a(t)dt as a unit of length and we 
take the energies as dimensionless (that is they are divided hj m < a >f /r^). 

We have considered: 

a) Harmonic oscillations 
That is to say: 

a{t) = I + 5 sin 27rt (2.3) 

This system resembles the so called full Fermi accelerator model. Such system was widely 
studied in simplified versions [0. Figure ^ shows the Poincare surface section for g = 
corresponding to the kinetic energy E versus the 6. Here we can distinguish three regions. 
The low energy region is fully chaotic, the intermediate energy region shows mixed dynamics 
(regular island surrounded by the chaotic sea) and finally the adiabatic region, where the 
dynamics become regular and we can see invariant curves. 

Now we assume that the moving wall is a finite barrier of potential of height Vq. So, the 
particles can leave the well and the limit between the bounded motion and the unbounded 
one is the line E = Vq. 

If Vq lies on the adiabatic region, that is there is an invariant curve below the limit line, 
the system will remain bounded. So to study the tdl we will set Vq less than the lower 
invariant curve. 

b) Sawtooth oscillations 



In this case: 

(l + 45t) if <t < 1/4 , 

«(i^) = \ (1 + 25) - 4(5t if 1/4 < t < 3/4 , (2.4) 

l-4(5)+4(5t if 3/4<t< 1 . 

The main feature of this system is that the periodic orbits are parabohc and non isolated. 
The two dimensional equivalent billiard is like billiards with neutral boundaries as polygonal 
ones !§,§]. For example, Fig. |^ shows a representative periodic orbit. Looking at the billiard 
as map T 

(g(0 = O),tana(0 = O)) ^ (g'(l), tana'(t = 1)) (2.5) 

the periodic orbit corresponds to a fixed point of T^. A straightforward calculation shows 
that T^ has eigenvalues Ai = A2 = 1 when it is linearized at this fixed point. So, it 
corresponds to a parabolic periodic orbit. On the other hand, given a fixed point (go, tanao), 
it can be shown that there exists an interval of g, / = {qmini Qmax), such that (g, tanao) is a 
fixed point if qel. That is, the fixed points are non isolated. 

III. THE DECAY OF CLASSICAL SYSTEMS 

To study the tdl, we fill a portion of the bounded region with A^o = 10^ particles whose 
initial conditions are uniformly (random) distributed in g and v = tana. We compute the 
ratio N{t)/NQ of the remaining population as function of time. 

Let us remark that the initial population must be chosen with low energy (in the chaotic 
region) to ensure the decay. If not so, particles can have initial conditions such that they will 
remain trapped because, in the case of harmonic oscillations, they can be on regular island 
included in the bounded region. On the other hand, we want to exclude the population of 
the region corresponding to very low energy to minimize the effect of the asymptotic orbits 
to the (trivial) parabolic subset of periodic orbits characterized hj E = 0. These trajectories 
are parallel lines to the time axis in the two dimensional pipe version of the system and they 



correspond to particles that having v = 0, they will never hit the boundary of the well (see 
appendix). So we have considered 

Figure ^ shows a part of the Poincare surface section of Fig. |^ and the region that was 
initially populated. We have taken Emin = 2.64 and Emax = 5.35. 

Figures |^ (a) and (b) show the results corresponding to the harmonic and sawtooth 
oscillations in log-lin and log-log plot respectively. We can see that the curves are quite 
different. In the case of Fig. |^(a) the log-lin plot shows an almost exponential behavior 
for the harmonic oscillations while Fig. ^(b) the log-log plot shows an algebraic one for 
the sawtooth oscillations. Such a difference can be explain in terms of the features of 
the periodic orbits that are included in the trapped region 0. In the case of harmonic 
oscillations, the trapped periodic orbits are hyperbolic and isolated (aside the subset E = 
that was not initially populated) while they are parabolic and non-isolated (like the bouncing 
ball orbits of the Sinai billiard) for the sawtooth oscillations. We stress that the deviation 
of the exponential behavior in the case of harmonic oscillations can be attributed to the 
population of very low energy region by indirect way (after collisions, see appendix). 

IV. QUANTUM SYSTEMS 

To study the quantum analogues of the systems, we want to solve the time dependent 
Schroedinger equation. 

^^^ = i/(g,t)|a,to,t> (4.1) 

being \a, to, to > an initial condition and H{q, t) the Hamiltonian operator H = T{p)+V{q, t) 
where: 

if < g < a(t) , 
V{q,t) = { -'- (4.2) 

Vo if a{t) < q . 



We have employed an expectral method (fast Fourier transform, FFT) based on the spht 
of the evolution operator U{t, to) |10[ defined by 



\a,t>=U{t,to)\a,to> (4.3) 



This technique was widely applied to solve bounded systems [11,12|, and transmission 



phenomena through moving barriers |13 



Such procedure does not avail with infinite barriers of potential, so to avoid the infinite 
barrier at g = 0, we have employed the symmetrized potential 

if < Igl < a(t) , 
V{q,t) = { I I- v; ^^_^^ 

[Vo ifa(t)<|g| . 

instead of potential ( [4.2| ) . 

We want to compute 

/■»(*) 
Pa(t) = \ <a,to,t\q> \^dq (4.5) 

J-a{t) 

that is the probability to detect the particle inside the well at time t. 

Because the system can be unbounded we will have a non stationary leaving flux of prob- 
ability going to g = ±oo. The spectral method automatically imposes periodic boundary 
conditions, so the flux of probability leaving at g = +{—)qmax will appear as coming at 
q = —{+)qmax- So we need to kill the leaving flux. To get this we have employed absorbing 



boundary conditions [M adding an absorptive (purely imaginary) static potential barrier 



y ^ JUO ^ iUo .^ gx 

cosh^ (7(g - qmax)) cosh^ (7(g + qmax)) 



to the actual (real) time dependent potential well |4.4| . Here, Uq and 7 are real parameters 
whose values were chosen to minimize the reflection and transmission coefficients of Va- All 
computations were performed taking a time step At = 0.00064 on a grid of 4096 steps of 
spatial sampling Ag = 0.00125. The absorbing region at the boundaries includes 250 spatial 
steps. This fact determines the value of 7. We have taken Uo = 14.79. Figures |^ summarizes 
the characteristic of the potential for the present calculation. 



On the other hand we have change the sharp step function in the potential |4.4| by a soft 
version hke a Saxon- Woods profile 

1 + expo(a(t)^ — g^j 



as was prescribed in Ref. [|T^ to improve the convergence of the numerical method. We have 
found that to ensure convergence and stability, we need that the potential vary from Vq/IO 
to 9Vo/10 in two spatial steps Ag. This determines the parameter b. 

V. THE DECAY OF QUANTUM SYSTEMS 

To populate the quantum system in an equivalent way as in the classical calculation, we 
consider the evolution of a quantum ensemble given by 

p(to) = ir^Zl l"'^o,io >< tt,io,io| (5.1) 

where N is the number of eigenstates of the static well of deep Vq = V{q,t = 0) whose 
eigenenergies lie on the region that we have populated the classical system. Such number 
depends on the value of h. Setting an appropriate value we calculate 

rait) ^ 1 

We have compute ( p.2|) taking h = 0.0255. For this value there are ten bounded eigenstates 
of the static well whose eigenenergies are such that Emm ^ E < E^ax- 

The results are shown in Figure ^ a) and b). The first shows a log-lin plot where we 
can see that both curves, corresponding to harmonic oscillations (curve A) and sawtooth 
oscillations (curve B) have non-exponential behavior. The second one shows a log-log plot 
that evidences an algebraic law from t > 25 for both curves unlike to the classical behaviors 
where for the same time interval they are qualitatively different (see Fig. ^(a) and (b)). 
Let us remark that even algebraic (oc l/t^), the characteristic exponents of the qtdl are 
quite different. The best fit for the qtdl corresponding to the harmonic oscillation system is 
f3 = 1.45 while it is /? = 1.11 for the sawtooth oscillation system. 
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P{t) ^ l< q\m\q >dq = -Y. Pc.it) . (5.2) 

J -at) JM tzii 



VI. SUMMARY AND CONCLUSIONS 

We have studied the time decay laws of two quasibounded systems in the framework 
of classical and quantum mechanics. The classical results show that the tdl's of the two 
systems are qualitatively different. One of the tdl results exponential while the other is 
algebraic. This dissimilitude can be explain in term of the diverse characteristics of the 
trapped periodic orbits in the systems. On the other hand, the quantum systems do not 
show such difference. We have obtained algebraic qtdl for both systems. However, the 
characteristic exponents of the qtdl are different. 

There are evidences that open quantum systems whose classical counterparts are chaotic 
while they remain bounded (with tdl almost exponential) can originate algebraic qtdl |l^ 



with different characteristic exponents. This effect can be explained by the distribution 
(Gaussian) of resonances, assuming an exponential decay of each resonance and an initial 
population of them. Varying the initial population, different characteristic exponents can 
be obtained. Such algebraic decay does not seem to be related to the existence of parabolic 
periodic orbits that generate algebraic tdl in classical systems [§]. In our systems, the initial 
population is the same for both systems so we think that the different exponents can be 
originate by the distribution of antiresonances (that is particular states that are refractory 
to absorb energy from the oscillating wall) into each well. 



On the other hand, in Ref. [Il6| , p!7| it is shown that, in general, the qtdl corresponding to 
chaotic systems have a power long time tail (oc 1/t^). The specific characteristic exponent f3, 
depends on the number of open channels for the decay. In particular, for one open channel 
the exponent is /3 = 1.5 [0. This fact, could be related to our result /5 = 1.45 for the 
harmonic oscillating wall system. Moreover, the characteristic exponent (3, corresponding 
to non-chaotic systems for one open channel and strong coupling between quasibounded 
states and the continuum is closer to /? = 1 [jl9[ because the widening of the distribution 



of resonance widths respect to those corresponding to chaotic systems. This result could 
be related to our result /3 = 1.11 in the case of sawtooth oscillating wall system. These 
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apparent connections are been studied. 

To finish, in spite of the classical systems and the quantum systems considered in the 
present work are not completely analogues (essentially because the smoothing of the quan- 
tum potential), we expect that the tdl corresponding to the classical systems with the 
smoothing potential will be the same as the step potential systems. As we remark in Sec. 
m, the tdl for the classical systems is determined by the linear stability of the trapped 
periodic orbits and such characteristic seems to be related to the time dependence (instead 
the spatial dependence) of the potential. 
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APPENDIX A: 

The present appendix is devoted to show the effect in the tdl when the very low energy 
region is initially populated. 

Let us consider a well of deep Vq. We have calculate the tdl starting off the uniform 



population of whole trapped region (that is < g < 1 and — a/SVq < v < \PWq. Fig. ||(a) 
shows the tdl as a function of time. Here, we can see an algebraic behavior, that is 

N{t)lN^ (X 1/t^ , (Al) 

where /? = 1. 

According Ref . , we expect a subset of periodic parabolic orbits included in the bounded 
region. This subset corresponds to particles that having f ^ generate trajectories which 
are asymptotic to the horizontal paths f = in the pipe version of our system. 

To show this, let us consider n(t) the fraction of initial conditions that spend a time t' 
greater than a given time t, to reach the oscillating wall. Looking at Fig. ||, we estimate, 
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for t > IJ^JWq, 



n{t) 



d{tanq 

- V -■ ^ -tan <t>min = i^+q)/t 



(A2) 



So that 



n{t) 



1 1 

X 



Following Ref. 



VWo t 
we assume for the long time tail 

pt'=oo 

N{t)/No ^ io / g{t')dt' 

Jt'=t 



(A3) 



(A4) 



where g{t)dt is the fraction of initial conditions for the first collision with the moving wall 
occurs between t and t + dt and u is the probability to leave the wall after one collision. As 

—dn 



9{t) 



dt 



(A5) 



we have the law lATI for the tdl. We stress that in the case of Fig. ^T] the algebraic tail 



becomes for t ~ 2/-\/2Vo but in the general, if we populate the low energy region such 
that |i;| < |fmaa;| < -\/2Vo wc cxpcct that the tdl exhibits a crossover between a stretched 
exponential and the algebraic decay law ^1| for long times as it shows in Fig. ^(b). 

When the initial population excludes the very low energy region, that is |fmm| ^ bl < 
kmaxi; wc havc uumerical evidence that the algebraic tail for long times follows 1/t^ (See 
Fig. |10D. So, it differs by one respect to the precedent. This difference also occurs in 
the decay of other quasibounded hamiltonian systems when a parabolic region is initially 
populated or it is populated by an indirect way in [^pl . 
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FIGURES 
FIG. 1. a) One trajectory in the two dimensional pipe version of the system. The coordinates 

are the time t and the position of the particle q. b) The same trajectory in the billiard version. This 

point of view is obtained by exploting the time lattice translation symmetry of the pipe setting 

9 = t — [t/T]T as a coordinate. 

FIG. 2. E (kinetic energy) vs. Q (reduced time) at q = 0. Poincare surface section correspond- 
ing to the system when the wall at q = a{t) oscillates armonically. Both, E and Q are dimensionless 
and the amplitude of the oscillation is 5 = 0.2 (in units of < a >t, see the text). 

FIG. 3. One parabolic periodic orbit corresponding to the sawtooth oscillation. Here q is the 
position of the particle while 6 is the reduced time. 

FIG. 4. E (kinetic energy) vs. Q at q = 0. This Poincare surface section shows only the 
low energy region and it was generated by one initial condition. The solid straight lines limit the 
initial population that is mentioned in the text. The dashed line indicates the deep of the well (it 
separates the bounded region and the unbounded region). E and are dimensionless. 

FIG. 5. a) Log-lin plot of the remaining population N/Nq inside the well against t (tdl). We 
have fixed r (the period of the moving wall) as the unit of time and < a >t (the mean value of 
the position of the moving wall) as the unit of length. (A) correspond to the harmonic oscillations 
system and it has an almost exponential behavior while (B) correspond to the sawtooth oscillations. 
In both cases the amplitude of the oscillation is 5 = 0.2 and the deep of the well is Vq = 8.63. b) 
The same curves of a) but in a log- log plot to show the algebraic behavior of (B). 

FIG. 6. This Figure shows the characteristic of the potential for the computation of qtdl. The 
spatial grid has 4096 steps. The absorbing region is marked by the shadow rectangles at the 
boundaries. Each side has 250 steps. 
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FIG. 7. a) Log-lin plot corresponding to the qtdl's P vs. t. Curve (A) corresponds to harmonic 
oscillation. Curve (B) corresponds to sawtooth oscillations, b) Log-log plot for the same laws for 
t > 25 to show the algebraic behavior. The units are the same as Figure ^. 

FIG. 8. Graph to show the fraction of initial conditions for which the first collision on the 
moving wall occurs in times greater than t. Given an initial position q, the velocities are v = tanc/) 

such -(j)min < 4> < <t>max- 

FIG. 9. a) Log-log plot of the tdl corresponding to harmonic oscillating well system starting 
off the uniform population of whole trapped region. The deep of the well is Vq = 0.78 and 5 = 0.2. 
We have also drawn the straight line of slope one to clarify the characteristic exponent of the 
algebraic behavior. b)Log-log plot of the tdl corresponding to harmonic oscillating wall. In this 
case Vq = 19.7, 5 = 0.2 and the initial population is the same of a), that is \v\ < 1.25. Here we can 
see the algebraic behavior for the long time tail. We also draw the straight line of slope one. The 
unit of time and the unit of length are the same as in Figure ^. 

FIG. 10. Log-log plot of the tdl corresponding to harmonic oscillating wall. In this case 
Vq = 8.63, 5 = 0.2 and the initial population was uniform in the region shown in the Figure ^ 
which excludes the parabolic region. Here we can see the algebraic behavior for the long time tail. 
We also draw the straight line of slope 2. The unit of time and the unit of length are the same as 
in Figure |5[ 
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